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Abstract 

We study the neutraUno radiative decay into the lightest supersymmetric particle 
(LSP) in the framework of a minimal anomalous C/(l)' extension of the MSSM. It 
turns out that in a suitable decoupling limit the axino, which is present in the 
Stiickelberg multiplet, is the LSP. We compute the branching ratio (BR) for the 
decay of a neutralino into an axino and a photon. We find that in a wide region of 
the parameter space, the BR is higher than 93% in contrast with the typical value 
(< 1%) in the CMSSM. 
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1 Introduction 



A great deal of work has been done recently to embed the standard model of particle 
physics (SM) into a brane construction [1, 2, 3, 4]. This research is part of the effort, 
initiated in [5], to build a fully realistic four dimensional vacuum out of string theory. 
While the original models were formulated in the framework of the heterotic string, the 
most recent efforts were formulated for type II strings in order to take advantage of the 
recent work on moduli stabilization using fluxes. Such brane constructions naturally lead 
to extra anomalous [/(l)'s in the four dimensional low energy theory and, in turn, to 
the presence of possible heavy Z' particles in the spectrum. These particles should be 
among the early findings of LHC and besides for the above cited models they are also 
a prediction of many other theoretical models of the unification of forces (see [6] for a 
recent review). It is then of some interest to know if these Z' particles contribute to the 
cancellation of the gauge anomaly in the way predicted from string theory or not. In 
[7] some of the present authors have studied a supersymmetric (SUSY) extension of the 
minimal supersymmetric standard model (MSSM) in which the anomaly is cancelled d 
la Green-Schwarz. The model is only string-inspired and is not the low-energy sector of 
some brane construction. The reason of this choice rests in our curiosity to explore the 
phenomenology of these models keeping a high degree of flexibility, while avoiding the 
intricacies and uncertainties connected with a string theory construction. For previous 
work along these lines we refer to [8, 9, 10, 11, 12, 13, 14, 15, 16, 17]. In this work 
we continue the analysis of the axino interactions [18] , studying the neutralino radiative 
decay into the axino. The next to lightest supersymmetric particle (NLSP) radiative 
decay might be the first process where the LSP could be observed at LHC [19, 20, 21, 22]. 
In our model we assume the axino as the LSP and a generic MSSM neutrahno as the 
NLSP. Due to the axino interaction vertices the NLSP neutrahno can only decay into an 
axino plus either a photon or a Zq or SM fermions. We compute all the amplitudes for 
that decays and the BR for the decay of a neutralino into an axino and a photon. This 
is the plan of the paper: in Section 2 we describe our model. In Section 3 and 4 we find 
the LSP and study the axino interactions. In Section 5 we compute all the decay rates 
under study and finally in Section 6 we show the results and the related plots. Section 7 
is devoted to the conclusions. 



2 Model Setup 

In this section we briefly discuss our theoretical framework. We assume an extension of 
the MSSM with an additional abehan vector multiplet V^^^ with arbitrary charges. The 
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anomalies arc cancelled with the Green-Schwarz (GS) mechanism and with the General- 
ized Chern-Simons (GCS) terms. All the details can be found in [7]. All the MSSM fields 
are charged under the additional vector multiplet V^^\ with charges which are given in 
Table 1, where Qi, Li arc the left handed quarks and leptons respectively while f/f , _D?, Ef 
are the right handed up and down quarks and the electrically charged leptons. The su- 
perscript c stands for charge conjugation. The index i = 1, 2, 3 denotes the three different 
families. Hu^^ are the two Higgs scalars. 
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Table 1: Charge assignment. 



The key feature of this model is the mechanism of anomaly cancellation. As it is well 
known, the MSSM is anomaly free. In our MSSM extension all the anomalies that involve 
only the SU{3), SU{2) and U{1)y factors vanish identically. However, triangles with 
U{\y in the external legs in general are potentially anomalous. These anomalies are^ 

(1) 
(2) 
(3) 
(4) 
(5) 

All the remaining anomalies that involve f/(l)'s vanish identically due to group theoretical 
arguments (sec Chapter 22 of [23]). Consistency of the model is achieved by the contribu- 
tion of a Stiickclbcrg field S and its appropriate couplings to the anomalous U{1)' . The 
Stiickelberg lagrangian written in terms of superfields is [24] 
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^We are working in an effective field theory framework and we ignore throughout the paper all the 
gravitational effects. In particular, we do not consider the gravitational anomalies which, however, could 
be canceled by the Green-Schwarz mechanism. Moreover it can be shown that A^^^ =0 [7] . 
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where the index a = 0, ... ,2 runs over the U{iy, U{1)y and SU{2) gauge groups respec- 
tively. The Stiickelberg multiplet is a chiral superfield 

n - 1 - 

and transforms under U{1)' as 

y{o) ^ 1/(0) +^ (A _ At) 

S 5 - 4z 63 A (8) 

where 63 is a constant related to the Z' mass. In our model there are two mechanisms 
that give mass to the gauge bosons: (i) the Stiickelberg mechanism and (ii) the Higgs 
mechanism. In the following we assume^ that 



Qh^ = Qh, = 

The mass terms for the gauge fields are given by 



(9) 



1/(1)/^ 



(10) 



with being the gauge boson mass matrix 
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9 ^,2 „.2 



(11) 



92 



2v^ 



where My(o) = Ab^go is the mass parameter for the anomalous U{1) and it is assumed to 
be in the TeV range. The lower dots denote the obvious terms under symmetrization. 
After diagonalization, we obtain the eigenstates 



9l 



and the corresponding masses 
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(13) 
(14) 

(15) 
(16) 
(17) 



''We impose this condition to simplify our computations and to give a compact analytical expressions. 
There are no obstructions to set Qh^ ^ 0, Qh^ 7^ 0. 
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Finally the rotation matrix from the hypercharge to the photon basis is 
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V ^ J 



y/gf+gl y/gf+gl I 
where i, j = 0, 1, 2. 

We now give the expansion of the lagrangian piece defined in (6) in component 
fields only for the part that is needed in the following sections. Using the Wess-Zumino 
gauge we get 

2 

r 



^ a=0 



"2 



AWaW^°V5 + (0^1) 



+ h.c. 



(19) 



As it was pointed out in [8], the Stiickelberg mechanism is not enough to cancel all the 
anomalies. Mixed anomalies between anomalous and non-anomalous factors require an 
additional mechanism to ensure consistency of the model: non gauge invariant GCS terms 
must be added. In our case, the GCS terms have the form [25] 



+ 



+dGTr 



(l^(2)i5-y(o) _ y(o)D«^(2)) p^W + n.a.c + h.c 



(20) 



6*252 



where n.a.c. refers to non abclian completion terms. The b constants in (6) and the d 
constants in (20) arc fixed by the anomaly cancellation procedure (for details see [7]). 

For a symmetric distribution of the anomaly, we have 

M'^ ,rn, ^« .(2). A^'^ 
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(21) 



It is worth noting that the GCS coefficients d4_5^6 are fully determined in terms of the ^'s 
by the gauge invariance, while the b^2^^^ depend only on the free parameter 63, which is 
related to the mass of the anomalous U{1). 
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The soft breaking sector of the model is given by 

>C.o/. = Cr^^r - I (MoA('^)aW + h.c.) - I (^^i^si^s + h.c)j (22) 

where C^^^'^ is the usual soft susy breaking lagrangian while A*^°^ is the gaugino of the 
added U{iy and ips is the axino. The axino soft mass term deserves some comment: 
from [26] we know that a fermionic mass term for a chiral multiplet is not allowed in pres- 
ence of Yukawa interactions in which this chiral multiplet is involved. But in the classical 
Lagrangian the Stiickelberg multiplet cannot contribute to supcrpotential terms given that 
the gauge invariance given from our U(l)' symmetry (8) requires non-holomorphicity in 
the chiral fields. In fact in our model both the axino and the axion couple only through 
GS interactions. It is worth noting that a mass term for the axion is instead not allowed 
since it transforms non trivially under the anomalous 11(1)' gauge transformation (8). At 
first sight our lagrangian may look not the most general possible one. In particular, an 
explicit Fayet-Iliopoulos (FI) term ^V^^^ and an explicit kinetic mixing term 
could be added. For what concerns the FI term, it is well known that in certain string- 
inspired models (see, e.g. [27, 28]), a one loop FI term is absent, even if Tr{Q) ^ 0. 
This is in apparent conflict with the observation [29] that in field theory a quadratically 
divergent FI term is always generated at one loop. The solution to this paradox is that 
in the low-energy lagrangian there should be a counterterm, which compensates precisely, 
i.e. both the divergent and the finite part of the one loop contribution. We do not write 
explicitly this counterterm, since its exact expression is model and regularization depen- 
dent, but we implicitly assume that such a cancellation occurs. For what concerns the 
kinetic mixing term, also this arises at one loop level and 8 oc Tr((5y). To simplify our 
computations we cancel the b term by choosing 

Qq = Ql (23) 

Moreover the constraint (9) implies A^^'' ~ 0, cancelling another possible source of ki- 
netic mixing. Finally, as mentioned before, we do not take into account the anomaly 
cancellation in the gravitational sector. 

3 Neutralino Sector 

Assuming the conservation of R-parity the LSP is a good weak interacting massive particle 
(WIMP) dark matter candidate. As in the MSSM the LSP is given by a linear combination 
of fields in the neutralino sector. The general form of the neutralino mass matrix is given 
in [7]. Written in the interaction eigenstate basis (■0°)^ = {i^s, A*^^\ Xf\ h!^, hi) it 
is a six-by-six matrix. From the point of view of the strength of the interactions the two 
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extra states are not on the same footing with respect to the standard ones. The axino and 
the extra gaugino A*^'^-' dubbed primeino are in fact extremely weak interacting massive 
particle (XWIMP). Thus we arc interested in situations in which the extremely weak 
sector is decoupled from the standard one and the LSP belongs to this sector. This can 
be achieved at tree level with the choice (9). The neutralino mass matrix becomes 
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(24) 



where Mg, Mq, Mi, M2 are the soft masses coming from the soft breaking terms (22) 
while My(o) is given in (11). It is worth noting that the D terms and kinetic mixing terms 
can be neglected in the tree-level computations of the eigenvalues and eigenstates. 

Moreover we make the assumption that Mq S> Ms, My(o) so the eigenstates are nearly 
pure axino and primeino, and we suppose that the axino V's is the LSP [18]. We consider 
the NLSP to be either a mixture of the bino A*^^^ and higgsino or a nearly pure wino Ag^^ 
The first situation is a typical configuration of the mSUGRA parameter space^ while 
the second situation is naturally realized in anomaly mediated supersymmetry breaking 
scenarios. 



4 Axino Interactions 

The axino interactions can be read off from the interaction lagrangian (19). The relevant 
term, written in terms of four components Majorana spinors^, is given by 

where the h^2^ coefficients are given in (21). The related interaction vertex is depicted in 
Fig. 1 and the corresponding Feynman rule is 

C(«)[7^7l^A:^ (26) 
where is the momentum of the outgoing vector and C^^^'s 

C« = ^glbf (27) 



5 Or in the so called Constrained MSSM (CMSSM) 
^The gamma matrices 7*^ are in the Weyl representation. 
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(a) 



Figure 1: Gaugino-axino- vector interaction vertex. 



are factors which contain the coupling constants and the parameters 63 which are related 
to the anomalous 11(1)' [7]. Therefore C^""^ <^ ga and the axino interactions will be 
extremely weak, being suppressed by an order of magnitude factor with respect to the 
weak interactions. 



In our model there are only few allowed decay channels. By assuming R parity conser- 
vation, the NLSP can decay only into the LSP plus other SM particles. Since the axino 
(LSP) interacts only with vertices of Fig. 1, only the gaugino fraction of the neutrahno 
(NLSP) gives a contribution. It can decay at the leading order only into gauge vector 
bosons or SM fermions (see Fig. 2). Moreover, since we assume a near mass degeneracy 
between ips and the NLSP, the production of the tt pair is suppressed, while the pro- 
duction of the Zq is allowed only for high neutralino masses. To emphasize the fact that 
only the gaugino component plays a key role in the radiative decay, from now on we will 
refer to gaugino instead of neutralino. For the cases in hand (bino-higgsino or pure wino 
NLSP) the BR of the process will be independent on the gaugino fraction since this is 
factorized in all the amplitude. 

5.1 A^"^ — )• IPs'! decay rate 

The corresponding Feynman diagram is depicted in Fig. 2a. The decay rate is given by 



5 Gaugino Decay Channels 




(28) 



where 




M!-Ml 

2Ma 



(29) 
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(a) 

Figure 2: Gaugino decay modes, (a) Axino- photon production (b) Axino-Zo production 
(c) Axino-// pair production. 

is the purely spatial momentum of the outgoing photon in the gaugino rest frame, Ma is 
the gaugino mass and Mg is the axino mass and 

M^:^ = ^C(«)^Z575 [Tm, 7.] u^k'^ (30) 

is the Feynman amplitude of the process with the outgoing photon momentum. The 
following coefficients 

C^^^ = ^9lh?^si^ew (31) 
are the factors C^"^ rotated on the photon eigenstate. The square amplitude is given by 



2 

A 



a 
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-V'"'k''k>^T^,,p (32) 
where 

T^uap = Tr [7m, 7.] J^„ b., 7p] ] - MsM^Tr [ [7^, 7,] [7,, 7,] ] (33) 
with ps and Pa the axino and gaugino momenta respectively. The pre-f actor 1/2 is the 
average over the spin states of the gaugino. Performing the computations we finally get 

r(«) = UA^'^^? (AM)3(AM + 2M,)3 

^ i^o^ j 5„^^ 10247r5(AM + M5)3M|, ' 

where AM = — Ms, Mz' the Z' mass, A^'^^ is the anomaly factor and 

R^^^ = cos' ew 

= sin' % (35) 
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5.2 A^'') ^ ^1)sZq decay rate 



The corresponding Feynman diagram is depicted in Fig. 2b. The situation is very similar 
to the previous case except that we have a massive vector in the final state. The decay 
rate is given by 

where 



m; -Ml + Ml 



is the purely spatial momentum of the outgoing Zq in the gaugino rest frame. The square 
modulus of the amplitude is given by 



{-V"'' + ^) k'^k'^T.^ap (38) 



where 



= -V2glb^^hmew 
eg = ^^2^4'^ cos (39) 



are the factor C^°'^ rotated on the Zq eigenstate while T^i^o-io is given in (33). The pre- factor 
1/2 is related to the average over the spin states of the gaugino. We finally get 

p(a) _ r 4 o(a) [ [i^Mf - M|J ((AM + 2M,)^ - Ml) ' 



X 



{2 (AMY + Ml) 
20 487r5( [{AMf - Ml] M|, 

where AM = Ma — Ms, Mz' the Z' mass, A^""^ is the anomaly factor and 

i?2=sin2% 

i?|; = cos2eH^ (41) 



5.3 A^'^^ '^sff decay rate 

The decay is kinematically allowed only if AM > 2m/. The corresponding Feynman 
diagram is depicted in Fig. 2c. The decay rate is given by 



max m£ix 



// (27r)-^32M3 7^.^.in ''Jmlr^ ff\ ^ ^ 
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where rriu and 777.23 are kinematic variables defined as 



m?2 - 



= Ml+m)-2MaEf (43) 
= Ml + Ml-2MaEs (44) 



and the extremes of integration are 



2 max ^ J\ /f ™ ^ 2 

12 
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(M„ - mfY (45) 
= (M5 + m/)' (46) 



m 



2 max 
23 



m: 



,2 min 
23 



where 



= (^2* + ^3)' - {^Ef -mj- ^Ef - mj) (47) 
= (£;,* + E*,)' - [^Ef -m} + ^Ef - mjf (48) 



E; = — ^ (49) 

E* = ^ ^ 50 

^ 2mi2 ^ ^ 

and m/ is the fermion mass. The Feynman amphtude of the process is ^ 

- i^-fC^zlk'-us^^ [7., 7.] «A^^3^«/7p(^f - af 75)^/ (51) 



Zq 



where gj, and are respectively the electric charge, vectorial and axial coupling of 
the Zq, C^^ and C^"^ are given respectively in eq. (31) and (39). More details can be 
found in Appendix A. 



6 Gaugino Radiative Decay. Branching Ratio 

In this section we summarize the results about the BR of the gaugino radiative decay. 
We consider an axino LSP and the NLSP to be either a bino-higgsino mixture or a nearly 
pure wino. We solve numerically the integral (42) and we express the result as a function 
of the axino mass Ms and the mass gap AM /Ms- The axino mass ranges from 25 GeV up 

''In the i?| gauge, ^ = 1. 
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Figure 3: Some diagrams contributing to NLSP decay in the MSSM: (a) NLSP LSP 7 
(b) NLSP ^ LSP // 

to 2 TeV while AM/Ms ranges from 0.5% up to 5%. The results are shown in the contour 
plots in Fig. 4-6 as functions of Ms and AM/Ms- There is no dependence on Mz' and 
on the anomaly involved in the process A''°'\ since they factorize in the decay rates in b^^ 
(see eq. (21)) and they cancel out in the BR computation. It is worth noting that there 
is no dependence on the gaugino fraction of the neutralino. As expected the BR is very 
high both for the bino-higgsino and wino case and the corresponding plots (Fig. 4) have 
no substantial differences. In a wide region of the parameters the BR is higher than 93% 
since the contribution coming from Fig. 2b is kinematically forbidden and the correction 
coming from Fig. 2c is only few percents (a second order process in perturbation theory). 
The situation is very different from the CMSSM case where the BR is lower than 1%, so 
it is never dominant [19]. In the unconstrained MSSM we expect the one loop process 
N2 Ni'-f (see fig. 3(a)) to be suppressed with respect to tree level process N2 -> Nif f 
(see fig. 3(b)) although the tree level decay contribution can be lowered by a suitable 
choice of the free parameters. 

The high AM regions plotted in Fig. 5 deserves a comment. In this region the contri- 
bution coming from Fig. 2b now is kinematically allowed but suppressed because AM is 
not very high with respect to M^^. In this case an important role is played by the rotation 
factors (35) and (41). In the bino-higgsino NLSP case they favor the radiative decay, while 
in the wino NLSP case they favor the decay with the Zq production and so there is and 
interplay between the kinematic suppression and the rotation factor enhancement. The 
result is that for the bino-higgsino NLSP case, the radiative process is still the dominant 
one, while for the wino NLSP case the decay rate for Zq production can be higher than 
the radiative decay rate. 

The WMAP allowed regions of parameters [18] are plotted in Fig. 6. The white region 
represents either BR< 93% or the WMAP forbidden region. Only for the wino case we 
have a WMAP allowed white region which is a tiny vertical strip at Ms ~ 1.6 TeV and 
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Figure 4: Branching Ratio for the gaugino radiative decay. On the x-axis we have the 
axino mass, while on the y-axis we have the percentage mass gap. 
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Figure 5: Branching Ratio for the gaugino radiative decay. Zoom in the high AM region. 
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Figure 6: Branching Ratio for the gaugino radiative decay. WMAP allowed region. 
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Wino Production 
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Figure 7: Number of wino produced inside the WMAP allowed region. Black scale means 
more than 10'^ then we lower the order of magnitude by a factor of 10 down to 1 (white 
scale) . 

AM /Ms > 0.04. We see that there is a huge difference in the allowed region for the 
two cases. The bino-higgsino case is allowed at low masses {Ms < 700 GeV) while the 
wino case is allowed at high masses {Ms > 700 GeV). In both cases the vertical strip 
Ms < 50 GeV is excluded by the lower mass limit on the MSSM neutralinos [30] and the 
Zq production region is forbidden. 

By assuming a center of mass energy of 14 TeV and an integrated luminosity of 
100 fb^^, we give an estimate of the number of NLSP produced. We focus on the wino case. 
The leading production processes are [31]: qq — ?■ A2A2, qg — )■ \2Q, QQ ^29, QQ' ^2^f, 
where Af is the wino-like chargino. The parton cross sections were integrated using the 
parton distribution function (PDF) package [32]. As an example we fixed the percentage 
mass gap at 5% and the gluino mass is M3 ~ 3M2 (typical in anomaly mediation). The 
wino-like chargino mass is the wino NLSP mass, since the mass degeneracy between the 
two states. As a simplifying assumption, we considered an universal squark mass niq for 
all the squarks involved. The dominant process is qg — )■ A2? because of the gluon PDF. 
The number of NLSP produced is given in Fig. 7 as a function of the wino mass M2 
and the squarks mass rriq. Since the BR is almost close to one this is also the number of 
photons in the final state. The number of events is always greater than 10 while in the low 
mass region the number can be greater than 1000. It should be stressed that our results 
are not a distinctive signature of a SUSY anomalous 17(1)' model but they are rather a 
consequence of the absence of a direct coupling of the LSP with (s)fermions. However 
this result can be combined with direct [7] and indirect searches [13] of anomalous Z' 
triangle interactions, measurements of the fermionic couplings by Z' decay width and 
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forward-backward asymmetries at LHC. We postpone a detailed analysis of these points 
in a forthcoming paper [33]. 

7 Conclusions 

We computed the BR for the radiative decay of the NLSP in the model described in [7]. 
Motivated by [18], we considered an axino LSP and the NLSP to be alternatively a bino- 
higgsino mixture or a nearly pure wino. In both cases wc found a very high BR (> 93%). 
This result is different from the CMSSM case where the BR is typically very low (< 1%). 
The corresponding WMAP allowed region are very different: low axino masses {Ms < 700 
GeV) for the bino-higgsino case, high masses {Ms > 700 GeV) for the wino case. This 
result could be used to discriminate between the two options for a possible decay in LHC. 
Anyway we will come back to this point in a forthcoming paper where we analyze in 
details the number of events for the gaugino radiative decay inside LHC [33] . 
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A A^") i^sff decay rate 

The amplitude in given by eq. (51). The corresponding square modulus is 
\Mfl\' = -32 

with Qf the electric charge of the fermion, e the corresponding coupling constant, Vf{af) 
the vectorial (axial) coupling of the fermion to Zq, qz^ the corresponding coupling con- 
stant, ruf the fermion mass, k the virtual gauge boson momentum, C'^/^o S^^®^ by (31), 
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(52) 
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(39) and 



% = m){paPs) - M,Ms 2m) + 3{pfPf)m} + {pfpj) 



+ 
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-'-a — 



-(P/P/) {PaPf){PfPs) + {PaPf){pfPs) +rn) ipaPs){PfPf) + 
-'^{PaPf){PfPs) - {PaPf){PfPs) - {PaPf){PfPs) " '^{PaP f){p fPs) 

{PaPf-){PfPs) + {PaPf){PfPs) mj + MiMs [mj - (p/P/)^] + 

-(P/P/) {PaPf){PfPs) + {PaPf-){PfPs) 



(53) 



where Pa, Ps, Pf and pj are the gaugino, axino and SM fermions 4- momenta respectively. 
The factor 32 contains also the 1/2 factor coming from the average of the gaugino spin 
states. The scalar products of the momenta are given by 



PaPS = 


MaEs 


PaPf = 


MaEf 


PaPf- = 


Ma{-Ef -Es + Ma) 


PfPS = 




PfPf = 


^{M^ + Ml-2EsMa-2m}) 
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\ {Ml -Ml- 2EfMa) 




Ml + Ml-2EsMa 



(54) 



and the integration variables are introduced by the equations 

M^ + m'j — m\2 



Ef — 
Es — 



2Ma 
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(55) 
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